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Abstract 

A study is made of the behavior of unstable states in simple models 
which nevertheless are realistic representations of situations occurring in 
nature. It is demonstrated that a non-exponential decay pattern will 
ultimately dominate decay due to a lower limit to the energy. The survival 
rate approaches zero faster than the inverse square of the time when the 
time goes to infinity. 

1 Introduction 

In this article a study is made of the nature of the decay of unstable states in 
a nonrelativistic setting. One would expect this decay always to be exponen- 
tial after some time, like the radioactive decay one finds in nature. It is the 
purpose of this article to demonstrate in a simple model with a ground state 
that next to a possible nondccaying bound state contribution there will be a 
non-exponential contribution to the decay pattern which will ultimately dom- 
inate the exponential decay. In practice this non-exponential contribution is 
extremely small and under normal conditions undetectable. Moreover it will 
arise only when one is able to keep a state coherent for a long time, which is 
not the case normally. This phenomenon is not unknown and is discussed in 
general terms in the literatureQ but a model demonstration may be helpful for 
its understanding. 

Starting point is a simplified model giving rise to exact exponential decay. This 
will be worked out in Section 2. In Section 3 the model is modified as to make 
it more realistic. This modified model nevertheless leads to a superposition of 
exponential decay contributions of which one will ultimately dominate. Both 
models have an essential flaw which is caused by the fact that no absolute lower 

^J.J. Sakurai, "Modern Quantum Mechanics, Revised Edition", Supplement II (Addison 
Wesley, 1994); G-C. Cho, H. Kasari and Y. Yamaguchi, Prog. Theor. Phys. 90, 803, (1993) 
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limit to the energy has been imposed. If an energy lower limit is taken into 
accomit then a non-exponential term appears to be unavoidable. Section 4 
gives a discussion of this effect. 



2 A simple, exactly solvable model 

Consider a free Hamiltonian Ho{= Hq), a set of kets and a particular ket 

I a) with the properties: 

Ho\E)=E\E), Ho\a)=a\a), -c» < < c» (1) 

{E\E') = -6(E- E'), (E\a) = (a\E) ^0, (a\a) = 1 (2) 

P 

/OC 
\E){E\dE+\a){a\=I, p>0 (3) 

-oo 

Here p is the density of free energy-eigenstates, assumed to be independent 
of E. Furthermore, there is an interaction Hamiltonian H'{= H''^) with the 
properties: 

/oo 
\E)dE (4) 
-oo 

Here we take a to be real, positive and independent of E. We have: 

{a\H'\E) ^ {E\H'\a) ^ a (5) 

This is the transition matrix element. 

Let the total Hamiltonian H be defined as: 

H = Ho + H' (6) 

We are interested in expressions of the form: 

{a\e-*"\a) (7) 

These are transition matrix elements for the transition of a state at time t — 
to the same state at time t > 0. Traditionally, one reasons as follows: Up to 
time t = the total Hamiltonian is Hq which allows an experimenter to prepare 
a system in a pure energy eigenstate \a). Then, at time t — the "perturbation" 
H' is "switched on" and at a later time t "switched off' . From that time on one 
can carry out an analysis of the resulting state, in particular one can try to find 
out what the chances are that the original state is found back. According to 
Heisenbergs uncertainty principle for time and energy it takes an infinite amount 
of time to prepare the original state and to analyse the final results, but the 
available time to do so is indeed unlimited. The method of time dependent 
perturbation theory leads then to the famous "Golden Rule" of Fermij^ 

^See any textbook on quantum mechanics. 
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In order to evaluate these matrix elements we put H in its spectral form: 

/oo 
A|A)(A|dA 
-OO 

H\\) = A|A) 
(A|A') = ^(A-A') 



\X){\\d\ = I 

Here we have made the assumption that H has no discrete eigenkets. 
We make the following expansion: 

/oo 
fx{E)\E)dE + cx\a) 
-co 

The normalization condition then gives: 

1 r°° 

- / mE)fx' {E)dE + clcy = S{X - A') 

P J —oo 

Next we solve the equation: 

{Ho + H' - X)\~X} = 

We have: 

{H„ + H'-X)\X) 

/oc 
fx{E)\E}dE + cx\a)] 
-OO 

/OO 
fxiE)[iE-X)\E)+a\a)]dE 
-OO 

/OC 
\E)dE] = 
-OO 

This can be satisfied only if: 

{E-X)fxiE)+pacx = 

and 

/OO 
fx{E)dE = 
-OO 

so that 

fxiE) = -pacxP^^ + f}d{E - A) 

and 

(a - X)cx 



r ME)dE=-^^ 

J -OO 0, 



Prom this we find: 

- pac, r P-l-dE + = _ (20) 

j-oo ^ - A a 

The integral is zero and we find: 

(a-A)c. 
a 

If we insert this we find the following expression for fx{E): 

ME) = -c,[paP-J— + ^SiE ^ A)] (22) 
£j — A a 

Except for a phase factor the coefficient cx can be calculated from the normal- 
ization condition. We have: 



r^{E)h'{E)dE 

/OO i \ 

^[paP^-^ + ^SiE-X)] 

■[p'^PEh:i + ^^SiE~X')]dE 



/°° 1 1 



a -A' /■°° „,„ „ 1 



+ / 5{E - X').paP- — -dE 

a - X 

^i^^X}(o^rsiE-X).SiE-X')dE] 

^ J-oo 

The first integral on the right is not zero. We have: 



(24) 



r 1 1 

/OO poo 
iE-X-ie)iE-X'-ie) ''' " 7^ " ^^^^A^ 

6{E - X') ^ . - 7r2^(A - A') 

£;-A-ie 

Here the first integral is zero as follows from contour integration. The next two 
terms give: 



We have: 

Both terms together therefore give: 

2tt'^5{X - A') (27) 

We therefore find: 

/_>^-''^''^ = ''^*(*-^'> w 

The following two integrals together give: 

+ / 5{E - X').paP-^dE (29) 

= p(a - X)P-^ + p{a - X')P—^ 



X-X> ' X'-X 

The last integral is: 

(a- A)(a- A') p {a-Xf 



[ 6{E- X).S{E - X')dE = S{X - A') (30) 

J-00 ^ 



Everything taken together: 

/oo 
me)fy{E)dE 
-00 



/OO -1 A 

[paP—- + ^^6iE-X)] 
-00 ii - A a 

■[p^PEh:i + ^^^iE-X')]dE 



Aa-Xf 

Thus we find: 



cW{[^^^^ + pW]5{X-X')-p} 



(A|A') = - r fl{E)h'{E)dE + clcx' 

P J-00 



clcy [^^-^ + pa^TT^jSiX - A') = 5{X - A') 



and so we can make the following choice: 



-\/pa 
a — A + inpa^ 



(31) 



(32) 



= (33) 
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so that: 



ME) = + ^6{E - A)] (34) 

Finally we find the following exact expression: 

= ^ ^""i nn^ ^r P<^P-^\E)dE + ^|A) - \a)] (35) 
a — A + mpa'^ J-oo E — \ a 

The spectral representation of H becomes: 

/oo 
A|A)(A|rfA 
-OO 



.[pa / 

»/ — ( 



This immediately leads to the following expression: 



/C30 
-OO (a-A)2+7rVa^ 

/OO -1 \ 

/OO 1 \ 



Prom this it follows that: 



/OO i, - 

■^ (c,-A>'+.V«' ''^ 

This integral can be evaluated by contour integration. Wc have already assumed 
t > 0, therefore, by closing the contour with a semicircle at infinity in the lower 
half plane we find: 



f 



{a - A)2 + TrV^a* 

po^e^TT ^(A — a + iTTpa^) 



= — 27ri lim 



x^a-ivpa^ (A — a)^ + w'^p^a^ 



9 it \ 

= —Zm nm 



(39) 



Thus we find: 



x-a-iirpa^ X — a — inpa 
—2i'iTpa'' 

{a\e-^"\a) = g-f ("-»-''°') (40) 
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Figure 1: Formal survival rate as function of time. 

In the course of time the probability of finding the original state back is equal 
to: 

Wait) = \{a\e-^"\a)f = e"^* (41) 

This is the exponential law of radioactive decay. Note that the coefficient of t 
in the exponent is in agreement with Fermi's Golden Rule. 
Except for the simplified model specifications no approximation is made. Strik- 
ing is that for negative t the same integral over A generates a plus sign in the 
exponent, so that time symmetry is restored. There is apparently no question 
of time irreversibility. A simple kink in the time curve appears (see Figure 1). 

3 Extension of the model 

In Section 2 we assumed p and a to be independent of E. We now introduce 
E-dependence. We then have: 

-oo <E <oo (42) 
0, {a\a) = 1 (43) 

p > (44) 



Ho\E)=E\E), Ho\a)=a\a}, 
{E\E') = -^^5{E - E'), {E\a) = {a\E) 



r 



p{E)\E){E\dE+\a){a\=I, 
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and 

/oo 
p{E)a{E)\E)dE (45) 
-OO 

Again we define: 

/oo 
h{E)\E)dE + c^\a) (46) 
-OO 

and try to solve the eigenvalue equation: 

[H - A)|A) = {Ha + H'- A)|A) = (47) 
with the normalization condition: 

(A|A') - y_ -^^f*{E)fy{E)dE + clcy = 5{\ - A') (48) 
This now gives: 

/oo 
h{E)\E)dE + cx\a)] 
-OO 

' fx{E)[{E - X)\E) + a{E)\a)]dE (49) 

O 

/OO 
p{E)a{E)\E)dE] = 
-OG 

This leads to: 

{E-\)h{E)+p{E)a{E)cx = Q (50) 

and 

/OO 
a{E)fx{E)dE - (51) 
-OO 

So we find: 

A(£;) = ^p{E)a{E)c,,P^^ + /3<5(i? - A) (52) 



and 



CA / p(£;)a2(i;)P-^l^d£; + a(A)/3 = -(a- A)cA (53) 



Now the integral is not zero. We find instead what is called the Hilbert trans- 

form^of the function p{E)a'^{E). 

Define: 

7m^p{E)a\E) (54) 

It is assumed that rj{E) — + sufficiently fast when E ±oo. 
The Hilbert transform a{\) of rj{E) is defined as: 

- i r rj{E)P^^dE (55) 



•^Erdelyi, Bateman Manuscript Project, "Tables of Integral Transforms 11", page 243. (Mc- 
Graw Hill, 1954) 
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In terms of this function we then have: 

a — A — na{\) 
a(A) 

so we find: 

h{E) = -cME)a{E)P^^ + ^— - A)] 
In order to make use of the normaUzation condition we evaluate: 

^ -fmfx'iE) 



p{Ey 



The integral: 



11 11 

can be determined and gives (see the Appendix): 

The remaining terms give no problems and we end up with: 

(A|A') = -L.fi^E)h'{E)dE + clcy 

= ''^''^'^^f" + ^(A)}'^(A - A') = 5{\ - A') 

Apparently we can choose: 



C\ = 



a — \ — 7ri7(A) + mrj{\) 
This gives: 



|A) = ^ ^^^W- 7Tt[ r p{E)a(E)P^—\E}dE 

' ' a - A - 7rcr(A) + i7rr;(A) 7-00 E - ' 

a — A — -KoiX) , , , , 



Prom this it follows that: 



and thus we obtain: 

(a|e-#^|a) = / e-^^{a\X){X\a)d\ 



(65) 



[a- A-7rcr(A)]2 +7r2?72(A) 



This is again an exact expression. If a is small also a is small and a good 
approximation can be obtained by writing: 

_^//| f°° ri{a)e-'i^dX 
{a\e ^ \a) ^ ]_^^^_^_ ^^^^^p + ^2^2(^) (gg) 

= g-f [a-'r'^(a)-i7rT;(a)] 

Note that this is again in agreement with Fermi's Golden Rule. 
Let us now consider the exact integral expression. By making use of theorems 
on Hilbert transforms we can get insight into the analytical properties of the 
integrand. The relevant theorems are proven in the Appendix and sound: 

1. If <j{y) is the Hilbert transform of r]{x): 

1 r°° 1 

a(y) = - / jj(x)P dx (67) 

IT x-y 

then —r}{y) is the Hilbert transform of (j{x): 

1 f"^ 1 

n{y) = / a{x)P dx (68) 

7rJ_oo x-y 

2. The function: 

is zero in the lower half z-planc, analytic and regular in the upper half ^;-plane 
and has the property: 

lim £,{z) = c7(xo) + iT]{xQ) (70) 

Here x and .To arc points on the real axis. The function can be analytically 
continued from the upper to the lower half z-plane but may not be regular 
there. It is clear that ^* {£) has just the opposite properties. 

3. The function: 

1 /■°° 11 

aix,y) = - Tj{z)P P dz (71) 

T^J-oo z-x z-y 
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satisfies the property: 

a{x, y) = ^^^^ ~ ""^^^ + 7:rj{x)Six - y) (72) 
X y 

The latter theorem has been used before and will be used later on. 
We can now rewrite the matrix element in terms of ^, ^* and ry: 

/ , "H, ^ r 77(A)e-Tr^(iA 

<*''^ l°''y-. |.-A-<(A)l.|,-A-.e(A)] <'^' 

If r?(A) is regular in the lower half A-plane except for poles then since: 

and ^*(A) is regular in the lower half plane, also (the analytical continuation of) 
^(A) is regular in the lower half plane except for poles. For t ^ we can split 
the integral into two parts: 



„ [a-A-<(A)].[a-A-7re(A)] 
1 , r e-*^dA r e-T^dX 



(75) 



27ri J_^ a - A - 7r^(A) a- X-n^* (A) 

In order to evaluate this expression we have to solve the equation: 

a - A - <(A) = (76) 

It is interesting to see what happens when we introduce a scaling factor and 
replace r] by tt] (r > 0). Then ^ is replaced by r^. We then have to solve the 
equation: 

a - A - 7rr^(A) = (77) 

When r approaches zero one of the roots approaches a. Then in first approxi- 
mation we have: 

a - A - TTT^ia) = (78) 

and we find: 

Ao ~ a — TTT^^a) = a — T:T[a{a) + iri{a)] (79) 

This is what we have seen before. Since 77(a) > this zero lies in the lower half 
plane and since this root never becomes real for any value of r it stays in the 
lower half plane for all values of t. 

The other roots approach poles in ^(A). Let A' be such a pole. Then, near A', 
we have: 

and we find: 

Ao ~ A' + TTT-^ (81) 

a — A' 
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Since ^(A) is regular in the upper half plane, for sufficiently small r these roots 
are all located in the lower half plane and since the roots are never real for any 
value of T they remain in the lower half plane for all values of r. Meanwhile, 
^*(A) is regular in the lower half plane, so for sufficiently small r there cannot 
be solutions of the equation: 

a- A-7rrr(A) =0 (82) 

located in the lower half plane. Now when we let t move from infinitesimal to 
regular values zero's will not move in and cannot spontaneously be created so 
for any r there will be no solutions of the equation in the lower half plane. For 
t > the conclusion is that: 



t/ — ( 



a- A-7re(A) 
and we find: 



= (83) 



(a e-s-^ a) = — / r (84) 



Let Ao be a root. Then SjAq < and its contribution to the integral is: 

1 



7rC'(Ao) + 1 

We end up with the following expression for t > 0: 



(85) 



(a|e-*^|a) = ^ 7;,„e-*^^«e-*l^^°l (86) 

Ao 

where the 7 satisfy the necessary but not sufiicient condition: 

Ao 

The result is a sum of exponentials decreasing with time. One of them ultimately 
will dominate. Note that the time derivative of 

W„(i) = |(a|e-t^^|a)|2 (88) 

in the limit t | is not zero. Still there is time reversal symmetry. Therefore 
again a kink appears at f = and the time derivative at t = does not exist. 

4 The influence of the energy lower bound 

Without loss of generality we may assume that = is the lowest energy value 
of the free Hamiltonian Hq. In that case we have, with a assumed ^ 0: 

Hq\E) = E\E), iJo|a)=a|a), ^<E<oo (89) 
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{E\E') = ^^5{E-E'), {E\a) = {a\E)=Q, {a\a) = l (90) 

/•OO 

/ p{E)\E){E\dE+\a){a\=I, p>0 (91) 

Jo 

and 

/•OO 

H'\E) =a{E)\a), H'\a) = p{E)a{E)\E)dE (92) 

Jo 

As before we consider normalized eigenkets of the Hamiltonian H = Hq + H'. 
We put H in its spectral form: 

/OO 
A|A)(A|rfA (93) 
-OO 

F|A) = A|A) (94) 
(A|A') =5(A-A') (95) 

/OO 
|A)(A|rfA = / (96) 
-OO 

Note here that we do not assume a lower bound on the eigenvalues. We now 
write: 

nOO 

|A)= / fx{E)\E)dE + cx\a) (97) 
Jo 

The eigenvalue equation gives: 

{Ho + H'- A)|A) = {Ho + if' - A)[ H fx{E)\E)dE + cx\a)] 

Jo 

POO 

= / fx{E)[{E - X)\E) + a{E)\a)]dE (98) 
Jo 

/•OO 

+ cx[{a - X)\a) + / p{E)a{E)\E)dE] = 

from which it follows that: 

{E - X)fx{E) + p{E)a{E)cx = (99) 

and 

POO 

{a-\)cx+ / a{E)fx{E)dE = (100) 
Jo 

This gives, as before: 

fx{E) = -p{E)a{E)cxP^^ + I35{E - A) (101) 
When A > we obtain an equation for /3: 

-cxj^ v{E)P^^dE + a{X)p = -{a-X)cx (102) 
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which can be solved: 



^ = ^) ^'^^^ 

In that case we have: 

h{E) = -c,[p{E)aiE)P^ + '^-^^-^^^5{E - A)] (104) 

where: 

^^^^ "'^^^^E^'^^ ^ I JZ '?(^)^(^)^:^'^^ (105) 

and is therefore the Hilbert transform of the function ri{E)9{E). We define: 

fi{E) = v{E)e{E) (106) 
We get the following expression for the eigenkets of H: 

i^>=-«rsM*'"^'^^ (107) 

+ [a-E- Tra{E)]6{E - X)}\E)dE + cx\a) 

We evaluate: 

f ^/;(i=)/v(i=)<ii= 

= ^aCa' / -TJtM^^Pt^ + [a-E- na{E)]5{E - A)} (108) 

Jo mE) E-X 

.{ti{E)P-=^^ + [a-E- Tra{E)]S{E - X')}dE 
Here we have: 



dE 



= .'-^^ + .'v{x)S{x-x') 



(109) 



Again the remaining terms can immediately be evaluated and the result be- 
comes: 

= r -krJl{E)fx'{E)dE + c^cv 

= '^^'^^'^^KA) ~ ^ ~ ^ ^(A)}^(A - A') = -^(A - A') 
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for A, A' > 0. This allows the normalization constant to be determined except 
for a phase. We choose: 



^ -v^ ^ -v^ (111) 

a — A — 7rcr(A) + iwr]{X) a — A — 7r^*(A) 

Note the difference with the former expression. Moreover it is only valid for 
A > 0. 

Finally we obtain for positive A: 



a- A-7r^*(A) Jo «(-^) ^ - (^^2) 



a — A — Tra(X) , , , , , , 



and so we get: 



a — A — 7r^*(A) 

When A < it is not possible to solve for /3 because then 6{E — A) is always 
zero. In that case we have: 

A(E, = -£<|)5ffic, (114) 

When this is inserted we obtain: 

which gives: 

a - A - 7rCT(A) = (116) 

This is an equation for A to be solved. We have for A < that the integral and 
all its derivatives are positive so for any a with the property: 

a<lima-(A) (117) 

there is one and only one solution. If a is larger then there is no solution. It 
might be that this limit is positive infinite. Then there is always one solution. 
The solution is to be interpreted as a bound state. Let Aq be this solution. We 

have for A < 0: 
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(119) 



This gives: 

= CaCA' [tT + 1] 

From this it follows that: 

(Aol Ao) = cl^cxo [ira'iXo) + 1] = 1 (120) 

and thus we have: 

na'{Xo) 

where: 

1 v{E) 
nJo {E- Ao) 
The following expression results if there is a bound state: 



(a|e-^^|a)= / e-*^(a|A)(A|a)dA + e-^^«(a|Ao)(Ao|a) 
Jo 

77(A)e-^^dA e'TT^o 



(123) 



A-7rCT(A)][a- A-7ra*(A)] 7r(7'(Ao) + 1 



The second term on the right hand side oscillates forever. That means that 
when one prepares the state |a) this state may contain a contribution from a 
bound state which does not decay. The first term can be rewritten in the form: 

1 /-"^ e-^^dX r°° e-T^dX 

27ri7o a-A-7re(A) Jo a-X-'K^*{Xy ^ ' 

Since the integrands have a possible pole on the negative real axis, extension of 
the integration path to minus infinity is not immediately possible. We have to 
avoid this pole and we do that by choosing the path of integration for the first 
integral as shown in Figure 2. 

If the same integration path were chosen for the second integral the extension 
of the integration path to minus infinity would not have changed the results, 
because the integrands are equal. However, then the second integral would not 
be zero. This integral is only zero when the integration path be chosen as in 
Figure 3. 

The difference is just the contribution from the bound state pole and thus we 
find that: 
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X,-plane 


1 = 

/ 


bound state pole — 







Figure 2: Integration path of the first integral. 



?*,-plane 






;i = o 










bound state pole — ^ 





Figure 3: Integration path of the second integral. 
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with the path of integration taken along the entire real axis except for the bound 
state pole and which is passed through the upper half plane, i.e. the contour is 
that of Figure 2. The integrand is analytic and regular in the upper half plane 
(on this particular Rieniann sheet) except for the pole on the negative real axis. 
The right hand side of this expression causes trouble because the functions 
f]{X) and a{X) are definitely nonanalytic in the lower half plane and a contour 
integration in the way used above is not possible. Suppose that 77(A) be analytic 
and regular except for poles; r]{X) > and finite for A > 0; 77(0) finite. Then 
r]{X)9{X) is "analytic" in the sense that it has a branch cut which on both 
sides of the real axis stretches itself out towards infinity and passes through the 
origin. The cut separates two analytic functions, one of them being identically 
zero, the other being 77(A). By using a similar argument as before one can now 
easily prove that ^(A) is still regular in the upper half plane and that its analytic 
continuation is regular in the lower half plane except for poles and a branch cut 
starting at the origin. Now the point A = is necessarily a true branch point. 
We take the cut along the negative imaginary axis. Except for the bound state 
pole all poles are on the right hand side of the cut and in the lower half plane. 
In order to evaluate the integral we deform the contour as to wrap the cut: 
The right hand side is rotated clockwise and the left hand side anticlockwise 
towards the negative imaginary axis. In the mean time poles are passed which 
give exponential contributions to the integral. The problem now is to evaluate 
the remaining contour integral. In Figure 4 the branch cut, possible poles and 
the original integration path are drawn. Note how the branch point is evaded. 



X -plane 


integrati 


on path 


n 


9 


• *^ 

t branch cut 

• 





Figure 4: Intermediate integration path. 

After deformation of the contour the situation is as sketched in Figure 5. 
Let r be the final contour. The poles on the right hand side of the imaginary 
axis give damped oscillatory contributions, the pole on the left hand side gives 
a pure oscillatory contribution from the bound state. 
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X -plane 

^ r 


integration path F 




^ s 

/ 


branch cut 







Figure 5: Final integration path. 



The remaining contour integral is: 

■^^d\ 1 



2'Ki Jy ol ■ 



A - 7r^(A) 2m 
_ _ 1 



F{-iX)e-i^dX 



(126) 







where F{X) is the "jump across the cut" of the function: 



a - A - 7rC(A) 



(127) 



The integral is apparently a linear combination of an infinite number of decaying 
exponentials. The function F{X) cannot be identically equal to zero between 
and some point on the negative imaginary axis, otherwise A = is not a true 
branch point. We find therefore that for t ^ oo the integral will dominate any 
exponential of the type e~~ for r > 0. The conclusion is therefore that there 
will be a non-exponential contribution to the transition matrix element which 
will ultimately dominate any exponential decay. This is the anomaly. 
Let us now discuss some details of this anomaly. Let 77(A) be a real, rational 
analytic function of A, positive when A is positive. Then for negative X we find 
that: 

a{X) = --r,{X)\n{--) + -A{X,c) (128) 

TT C TT 

with c a positive number with the dimension of A. That branch is chosen where 
the logarithm is real. Then the function A{X, c) is a real and rational analytic 
function of A, finite around A = and which for negative A is sufficiently positive. 
Now we consider two cases: 
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1. ?7(0) > 0. Then for small A: 



a - A - 7r^(A) = a - A - na{X) « 77(A) H-^) (129) 



The "jump across the cut" of the integrand is then: 

1 1 

r?(A)ln(-|) ~r?(A)[ln(-A) + 27ri] 



- . A^ - ,,,n ^ X, . . . ^0 ^ (1^0) 



2. r]{0) = 0. We have: 

a — A — 7r5(A) = a — A — 7r(T(A) 



A (131) 

a-A + 77(A)ln( )-A(X,c) 

c 



and we find then for small A: 

F(A) = ^ 



a-A + r?(A)ln(-A)-A(A,c) 
1 



a - A + r?(A)[ln(-^) + 27rz] - A(A, c) 



(132) 



[a- A- A(A,c)]2 



(A-0) 



The conclusion is that the transition matrix element behaves non-exponentially 
but goes faster to zero than t^^ for large t. Correspondiugly the survival rate 
goes slower than exponential but faster than to zero for large t. 
Of course, the anomaly has been demonstrated only in a simplified model which 
however is realistic enough to warrant the expectation that the anomaly is 
characteristic for a much wider class of models. It may even be unavoidable. 
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A Appendix 

In this Appendix we prove the three theorems on Hilbert transforms mentioned 
in the text. 

1. Let a{y) be the Hilbert transform of ri{x): 

1 r°° 1 

- / r]{x)P dx = a{y) (133) 

J-oo x-y 

then —r]{y) is the Hilbert transform of a{x): 



(134) 



1 1 

- / a{x)P dx = -T]{y) 

IT J-oo x-y 

Proof: In the text we have already proven that 
1 Z""" 1 1 

^ / P P dz = 5{x - y) (135) 

TT^ J_oo x-z y-z 

We now have: 

rj{x)5{x — y)dx 

(136) 

= — / [ / Tl{x)P dx]P dz = v{y) 

J-oo J-oo x-z y-z 

and therefore: 

r]{y) = - a(z)P dz = / a(z)P dz (137) 

TTj-oo y-z ttJ-oo z-y 

2. The function: 

def 1 r a{x)+iri{x) 
^("^=2^y_ x-z ^'^^^ 

is zero in the lower half z-plane, analj^tic an regular in the upper half 
2;-plane and has the property: 

lim £,{z) = cf{xq) + ir}{xQ) (139) 
Here x and xq are points on the real axis. 

Proof: The function is obviously analytic and regular in the upper and 
lower half ^;-plane. Since the real axis acts as a closed branch cut the two 
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branches are not analytically connected. We have: 



zlxo 27ri zJ.a;o J _oo X — Z 

1 f°° a{x) + ijj{x) 



„ , . dx 

27ri J_oo X — Xo — le 

1 1 (140) 

= -r-^ / Wix) + iri{x)\.\P h 7ri(5(a; - xo)\dx 

2-Ki J_J- ^ ' ' X-Xq 

= ^[-»7(a;o) + i(^{x<i)] + ^[(t{xo) + iri{xo)\ 
= a{xQ) + i'n{xo) 

This proves the first part. We also have: 
linre(.) = ^lim r ^Ml^^x 



i\xo 27ri z}xo J_oo X — Z 

1 p + iri{x) 



2TTi 7-00 a; — + ie 



dx 



(141) 



1 /■°° 1 
= / [a(x) + ir]{x)].[P TviSlx — xo)]dx 

2TTt X- Xo 

= -^[-v{xo) + i<7{xo)] - ^[cr{xo) + iv{xo)] = 

This means that the regular function in the lower half plane is zero on the 
real axis and must therefore be zero throughout the lower half plane. 



3. The function: 



1 r°° 11 

a(x,y) = - ri(z)P P dz (142) 

T^J-oo z-x z-y 



satisfies the property: 



a{x,y) = ^^^^-^^ + TT7]{x)S{x-y) (143) 
x-y 



Proof: 



1 f°° 11 

- / zr](z)P P dz 

Ti-J-oo z-x z-y 



1 r°° 11 

- / \{z - x) + x]r](z)P P dz 

TTj-oo z-x z-y 

- / T](z)P dz + x- / r](z)P P dz 

ttJ_^ z-y J-oo z-x z-y 



(144) 



= a{y) + xa{x, y) = (j{x) + ya{x, y) 
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It follows from this that: 

{x - y)a{x, y) = a{x) - a{y) (145) 

and from this we obtain: 

a{x, y) = ""^""^ ~ ""^^^ + A{x)S{x - y) (146) 
x-y 

The A{x) can be evaluated: 

a{x) - a{y) 



/OO rC 
a{x,y)dy- / 
-OO J — ( 



x-y 



dy 



1 r°° 1 f°° 1 

-a{x) / P^—dy+ / <j{y)P— 



-dy 

y 



'Kr]{x) 



(147) 
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